AD-A168  885 
UNCLASSIFIED 


COMPARISON  OF  RAVLEIGH-TAVLOR  DROUTH  WITH  IS 
MICROMETERS  AND  825  MICROMETERS  LASERS(U)  NAVAL 
RESEARCH  LAB  HASHINGTON  DC  H  H  EMERV  ET  AL.  29  HAV  86 
NRL-HR-5796  F/G  18/1 


1/1 


Comparison  of  Rayleigh-Taylor  Growth 
with  1.0  jul m  and  0.25  /xm  Lasers 

M.  H.  Emery  and  J.  H.  Gardner 
Laboratory  for  Computational  Physics 
S.  E.  Bodner 
Plasma  Physics  Division 


This  work  was  supported  by  the  U  S,  Department  of  Energy  and  the  Office  of  Naval  Research. 


DTIC 

tUL'LECTE 

JUN  2  3  1986 1 


Approved  for  public  release,  distribution  unlimited 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b.  DECLASSIFICATION /  DOWNGRADING  SCHEDULE 


«  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NRL  Memorandum  Report  5796 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory 


6c  ADORE  SS  (City.  SWt,  ind  ZIP  Codt) 

Washington,  DC  20375-5000 


3  DISTRIBUTION! AVAILABILITY  OF  REPORT 

Approved  lor  public  release;  distribution  unlimited. 


5  MONITORING  ORGANIZATION  report  NUMBERCS) 


6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
(It  tpplictblt) 

Code  4040  Department  of  Energy 


7b  ADDRESS  (City.  Stitt,  ind  ZIP  Codt) 

Washington,  DC  20585 


Ba.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Department  of  Energy 


Sc  ADDRESS  (City,  Stitt,  ind  ZIP  Codt) 

Washington.  DC  20545 


8b  OFFICE  SYMBOL  19  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
Of  tppllttblt) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 

ELEMENT  NO  I  NO 


tl  TITLE  (Includt  Security  Qiaifkttion) 

Comparison  of  Rayleigh-Taylor  Growth  with  1.0  pm  and  0.25  pm  Lasers 


12  personal  author(S) 

Emery.  M.  11.,  Gardner,  J.  H.  and  Bodner,  S.  E. 


13a  TYPE  OF  REPORT  13b  TIME  COVERED  14  DATE  OF  REPORT  (Ytir.  Month.  Dty)  h S  PAGE  COUNT 

Interim  from  9  85  to  3/86  |  1986  May  29  I  14 


16  SUPPLEMENTARY  notation 

This  work  was  supported  by  the  U.S.  Department  of  Energy  and  the  Office  of  Naval  Research. 


17 

COSATI  CODES 

FIELD 

1  GROUP  ! 

1  SUB-GROUP 

18  SUBJECT  TERMS  (Continue  on  rtvtne  if  ntetsury  tnd  identify  by  block  number) 

Rayleigh-Taylor 

Numerical  simulation 

Laser-ablation 

Hydrodynamic  instability 

19  ABSTRACT  (Continue  on  raven*  it  necessary  and  Htntfy  by  block  number) 

-  It  is  shown  through  numerical  simulation,that  the  Rayleigh-Taylor  growth  rate  in  laser-ablatively  accelerated  targets 
is  reduced  below  the  classical  value  (kg)l.-,  by  factors  of  3-4  with  0.25  pm  laser  light.  The  simulation  results  are 
supported  by  an  analytical  expression  for  the  growth  rates.  These  results  provide  further  evidence  for  the  viability  of 
high-asneet  ratio  shells  in  direct-drive  laser  fusion. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  |21  A8STRACT  SECURITY  CLASSIFICATION 

□  UNCLASSIFIEO/UNUMITED  □  SAME  AS  RPT  □  OTIC  USERS  I  UNCL  ASSIFIED _ 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  |22b  TELEPHONE  (Inclodt  Arti  Codt)  22t  OFFICE  SYMBOL 

Mark  li.  Emery  1  (202 1  767-3196 _ Code 4040 _ 


OO  FORM  1473,  B«  MAR  83  APR  edition  may  be  uted  until  anhauited  SECURITY  CLASSIFICATION  OF  THIS  RAGE 

All  other  edition!  «r*  otMOiet* 


COMPARISON  OF  RAYLEIGH-TAYLOR  GROWTH  WITH  1.0  Mm  AND  0.25  nm  LASERS 


One  of  the  most  critical  issues  facing  direct  drive  laser  fusion  is  that  of  hydrody¬ 
namic  stability.  In  laser  fusion  a  pellet  is  imploded  by  ablating  material  from  the  outer 
surface.  The  acceleration  of  a  cold,  dense  shell  by  a  hot,  light  blowoff  plasma  is  an  un¬ 
stable  arrangement.  The  Rayleigh-Taylor  (RT)1,  or  interchange,  instability  transposes 
the  position  of  the  more  dense  fluid  with  the  fluid  of  lower  density  in  order  to  reach  a 
more  stable  configuration.  This  instability  is  a  potential  obstacle  to  laser  fusion  since 
it  can  cause  sheil  fracture,  fuel  mix  or,  in  its  mildest  form,  a  nonuniform  implosion  that 
will  reduce  the  convergence  and  severely  diminish  the  energy  gain  of  the  fusion  pellet. 

Perhaps  the  most  serious  implication  of  the  RT  instability  is  that  it  ultimately 
dictates  the  allowed  aspect  ratio  R/AR  (ratio  of  the  distance  pushed  to  the  shell 
thickness)  for  the  shell.  The  allowed  aspect  ratio  both  fixes  the  pressure  required  to 
drive  the  shell  inward  and  indirectly  determines  the  efficiency  of  the  energy  coupling 
between  the  laser  and  the  shell2.  If  the  RT  grows  at  its  classically  predicted  rate,  then 
thick,  low-aspect  ratio  shells  must  be  used.  This,  in  turn,  implies  laser  intensities  at  or 
above  the  threshold  for  plasma  instabilities  and  unacceptably  low  rocket  efficiencies. 


The  allowed  aspect  ratio  for  the  laser-ablative  case  is  a  function  of  the  RT  growth 
rate  and  the  wavelength  of  the  most  dangerous  mode.  Classical  theory1  predicts  that 
an  initial  perturbation,  tj0i  on  the  interface  of  an  unstably  stratified  fluid  in  planar 
geometry  will  grow  as  77(f)  =  rjae '7t,  where  7  =  (A  kg)1/2,  k  =  2t/A,  where  A  is 
the  wavelength  of  the  perturbation,  g  is  the  acceleration,  .-1  =  (ph  —  pi)  /  ( ph  -f  pi)  is 
the  Atwood  number  and  ph{pt)  is  the  density  of  the  heavy  (light)  fluid.  (A  ~  0(1) 
in  the  laser-ablative  case.)  If  the  perturbation  wavelength  is  much  smaller  than  the 
radius  (R0)  of  a  spherical  pellet,  then  the  planar  growth  rate  can  be  applied,  at  a 
first  approximation,  to  spherical  geometry.  Let  us  assume  that  the  growth  rate  for. 
the  laser-ablative  case  is  e(kg)1/2,  where  s  <  1.  If  the  shell  is  imploded  at  constant 
acceleration  half-way  inward,  a  good  approximation  to  present  designs,  then  the  number 
of  e-foldings  of  unstable  growth  is  7 1  =  r(2~i?c/A)1/2.  Classical,  incompressible  theory 
suggests  that  the  most  dangerous  mode  is  A  =  A R,  the  thickness  of  the  shell,  as  then* 
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the  inner  and  outer  surfaces  are  strongly  coupled.  However,  numerical  simulations  have 
shown  that  when  compressibility  and  nonlinear  effects  are  taken  into  account,  it  is 
those  perturbation  wavelengths  several,  0(3),  times  the  shell  thickness  which  are  most 
damaging  to  the  implosion  process3.  For  mathematical  convenience,  we  set  A  =  it AR, 
and  the  limiting  aspect  ratio  is  given  by  R0/ AR  =  l/2(^)2. 


This  limit  is  quite  severe  since  an  initial  aspect  ratio  of  10  or  so  can  become 
an  in-flight-aspect-ratio  (IFAR)  of  ~  150  as  a  result  of  compression  due  to  the  20-40 
Mb  ablation  pressure.  Consider  a  pellet  polished  so  that  after  compression  the  largest 
surface  perturbation  has  an  amplitude  of  —  100A  and  assume  that  the  implosion  process 
can  survive  6  e-foldings  of  RT  growth.  If  the  RT  grows  at  its  classical  rate  (e  =  1)  then 
the  maximum  IFAR  is  limited  to  a  value  less  than  20.  The  initial  aspect  ratio  for  this 
pellet  would  be  about  1.5  .  The  laser  intensity  required  to  generate  sufficient  pressure 
to  accelerate  this  shell  to  implosion  velocities  would  be  well  above  the  threshhold  for 
plasma  instabilities.  These  thick  low-aspect  ratio  shells  would  also  severely  limit  the 
efficiency  and  gain  of  a  fusion  pellet4. 


Both  numerical  simulations5  and  experiments6  at  l/ma  laser  light  indicate  that 
the  RT  growth  rate  in  the  la^er-ablative  case  is  about  1/2  of  the  classically  predicted 
value  (e  =  1/2).  This  would  allow  an  IFAR  ~  70,  which  is  marginal  at  best.  Unless 
the  RT  growth  rate  is  reduced  even  further  than  that  obtained  with  1/xm  laser  light, 
these  thin,  high-aspect  ratio  shells  will  rupture  well  before  the  fuel  compression  stage. 
It  is  this  enhanced  reduction  in  the  growth  rate  at  shorter  laser  wavelengths  we  wish 
to  discuss  in  this  letter. 


The  RT  instability  is  modeled  using  the  FAST2D  laser-shell  simulation  code.  This 
is  a  fully  two-dimensional  Cartesian  code  with  a  sliding  Eulerian  grid  with  variable 
grid  spacing.  The  grid  spacing  is  0.25/im  for  15  zones  on  either  side  of  the  ablation 
layer  and  increases  uniformly  to  2  /xm  for  most  of  the  rest  of  the  grid.  The  refined 
subzoning  follows  the  ablation  front  throughout  the  course  of  the  run.  FAST2D  solves 
the  ideal  hydrodynamic  equations  using  the  flux-corrected-transport  algorithms  with 
two-dimensional  classical  (T5/*)  plasma  thermal  conduction.  The  code  has  been  exten¬ 
sively  documented  against  experimental  data5,7  and  is  discussed  in  some  detail  in  Ref. 
7  and  references  therein.  The  initial  steady  state  density  profile  is  perturbed  at  its  peak 
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with  a  single  sinusoidal  mode  corresponding  to  a  total  initial  mass  perturbation  of  less 
than  0.5%.  The  laser  pulse  has  a  2  ns  Gaussian  rise  after  which  the  laser  intensity  is 

held  constant.  The  computational  growth  rates  are  obtained  by  Fourier  transforming 
the  summed  mass  of  the  foil  which  is  integrated  from  the  rear  of  the  foil  to  the  ablation 
edge  for  each  transverse  coordinate.  The  growth  rate  is  measured  during  the  constant 
phase  of  the  laser  pulse.  Typically  6-7  e-foldings  of  exponential  growth  are  obtained. 

Figure  1  compares  the  growth  rate  obtained  from  the  FAST2D  simulation 

( -  •  - )  to  (kg)1/2  ( - )  for  1/im  laser  light  with  an  absorbed  intensity  of 

1013  W l cm2  impinging  a  20/xm  thick  plastic  (CH)  target.  The  computational  growth 
rates  have  a  k1/2  dependence  with  a  moderately  strong  cutoff  for  the  short  wavelength 
modes.  The  computational  growth  rates  are  about  a  factor  of  1.5  -  2  less  than  the 
classically  predicted  value. 

The  inhibited  nature  of  the  RT  growth  rate  for  the  laser- ablative  case  is  even  more 
apparent  for  shorter  wavelength  lasers.  For  this  case,  the  laser  intensity  was  increased 
to  a  value  closer  to  fusion  conditions.  Figure  2  compares  the  FAST2D  simulation  results 

( - • - )  to  (kg)1/2  ( — )  for  l/4^m  laser  light  with  an  absorbed  laser  intensity 

of  3  x  1014W/c7n2.  The  target  is  a  20/um  thick  plastic  (CH)  foil.  The  numerical  growth 
rates  for  this  case  are  a  factor  of  3-4  less  than  the  classically  predicted  values.  The  ratio 
of  the  numerical  growth  rate  to  (kg)1/2  for  both  the  l  pm  results  (•,  from  Fig.  1)  and 
the  1/4 (xm  results  (x,  from  Fig.  2)  is  plotted  as  a  function  of  perturbation  wavelength 
in  Figure  3.  Note  that  for  wavelength  perturbations  of  O(60;um),  three  times  the  foil 
thickness,  the  growth  rate  has  been  reduced  by  a  factor  of  1.8  for  1/xm  lasers  and  by  a 
factor  of  3.5  for  l/4jum  lasers. 

The  strong  dependence  of  growth  rate  upon  the  laser  wavelength  can  be  understood 
in  terms  of  the  mechanism  that  we  believe  is  responsible  for  the  reduction  below  the 
classical  value  :  the  ablative  convection  of  the  vorticity  away  from  the  unstable  ablation 
surface8.  It  is  the  vorticity  generated  at  the  unstable  interface  which  controls  the 
interchange  of  the  two  fluids.  As  a  result  of  the  laser-ablative  process,  part  of  that 
vorticity  is  convected  away  thereby  reducing  the  growth  rate  of  the  instability.  If  the 
ablation  velocity  (Va)  is  included  in  the  two-dimensional  vorticity  equation  through  the 
convective  term,  the  growth  rate  for  the  RT  instability  in  a  laser  ablatively-accelerated 
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target  is  given  by9 


T  =  [(kVa/2)2  +  'To]1/2  -  kVaj2  (1) 

where  =  {kg)l/2 ,  the  classical  growth  rate.  When  the  ablation  velocity  is  large,  then 
the  classical  RT  growth  rate  is  strongly  reduced.  The  above  expression  (x)  is  compared 
to  the  FAST2D  simulation  results  (•)  in  Figures  1  and  2.  The  agreement  is  quite  good 
for  both  the  1/rm  and  the  1/4/im  cases.  The  ablation  velocities  are  taken  from  the 
numerical  simulations  and  are  measured  at  the  center  of  the  unstable  shear  layer  on 
the  ablation  front.  See  Figure  4.  The  ablation  velocities  have  the  values  of  1.3  x  106 
cm/s  for  the  1/xm  case,  at  1013W/em2,  and  5.5  x  10®  cm/s  for  the  1/4/xm  case,  at 
3  x  10 14W/cm2.  The  lack  of  agreement  at  the  short  wavelength  perturbations  for  the 
l^m  case  is  due  to  the  assumption  of  an  infinitely  thin  vortex  sheet  at  the  ablation 
layer.  If  in  Eq.  1  is  replaced  by  the  growth  rate  for  an  exponential  density  layer  of 
finite  thickness,  then  the  cutoff  is  reproduced9. 

The  strongly  inhibited  growth  for  l/4^m  lasers  is  due  to  a  very  large  ablation 
velocity.  Shorter  wavelength  lasers  deposit  their  energy  closer  to  the  ablation  front 
and  produce  a  higher  mass  ablation  rate.  This  higher  mass  ablation  rate  produces 
a  higher  acceleration  of  plasma  away  from  the  target  surface,  and  therefore  a  larger 
value  of  Va.  Figure  5  is  a  plot  of  the  ablation  velocity  (Va)  obtained  from  the  FAST2D 
simulation  model  as  a  function  of  the  distance  from  the  target  surface  for  three  laser 
wavelengths  (l,  1/2  and  1/4  jum).  All  the  calculations  were  with  a  20/mo  thick  CH 
foil  with  approximately  the  same  isentrope  with  an  absorbed  laser  intensity  of  5  x 
10 13W/ cm2.  Note  that  at  a  distance  of  f/xm  from  the  target  surface,  which  is  at  the 
center  of  the  unstable  shear  layer,  the  ablation  velocity  for  l/4/xm  light  is  over  three 
times  larger  than  for  l/im  light. 

This  result  can  be  seen  more  clearly  from  a  scaling  arguement.  The  expansion 
velocity  measured  near  critical  scales  as  V.  a  (7Aj?)1,3  and  the  ablation  pressure  scales 
as  Pa  a  (I ' \i)z'3  in  planar  geometry10.  I  is  the  laser  intensity  and  A*  is  the  laser  wave¬ 
length.  This  scaling  is  somewhat  different  for  spherical  geometry  (Pa  a.  (/  A^'  )2  3) 
11  but  the  argument  still  holds.  From  the  conservation  of  mass  paV’.,  =  i>zVc  where 


4 


pa(V4)  is  the  density  (velocity)  to  be  determined  very  near  the  ablation  layer.  Assume 
that  the  target  and  laser  profiles  can  be  designed  so  that  the  targets  are  on  the  same 
isentrope;  i.e.,  Pa/pa  =  constant,  irrespective  of  laser  intensity  and  wavelength.  If  the 
ablation  pressure  is  kept  constant  as  the  laser  wavelength  is  reduced  then  I  aX t  and 
the  ablation  velocity  scales  as  Va  a  A^1.  If  the  ablation  pressure  is  allowed  to  vary  then 
Va  a  AJl4//XS.  Thus,  although  the  blowoff  velocity  at  critical  decreases  with  decreasing 
laser  wavelength,  the  ablation  velocity  measured  near  the  target  surface  increases  with 
decreasing  laser  wavelength.  Note  that  the  target  acceleration  ( g  =  Pa/M,M  is  the 
target  mass)  is  constant  if  the  ablation  pressure  is  held  fixed  and  scales  as  g  a  A^2^3 
if  Pa  is  allowed  to  vary  with  laser  wavelength.  The  scaling  for  the  ablation  velocity 
is  a  stronger  function  of  laser  wavelength.  The  growth  rate  for  the  RT  instability  will 
then  be  more  strongly  inhibited  with  short  wavelength  lasers  than  with  l/xm  laser  light, 
since  the  vorticity  will  be  convected  away  at  a  faster  rate,  see  Eq.  1. 

We  have  shown  that  the  RT  growth  rate  for  laser  ablatively-accelerated  targets 
is  strongly  reduced  below  the  classically  predicted  value  for  short  wavelength  lasers. 
By  combining  the  most  damaging  perturbation  wavelength  of  O(60^m),  about  3  times 
the  shell  thickness,  with  a  growth  rate  that  is  3.5  times  smaller  than  the  classical 
value,  s  as  1/3.5,  and  up  to  6  e-foldings  of  exponential  growth,  we  conclude  that  the 
maximum  IFAR  for  a  reactor-sized  pellet  driven  with  l/'-ium  laser  light  is  P/AP  as 
220.  If  one  designs  for  a  direct  illumination  pellet  with  an  IFAR  of  order  150  using 
1/4/xm  light,  then  the  net  (rocket  times  absorption)  efficiency  can  be  as  high  as  10?o4. 
Pellet  gains  can  then  be  on  the  order  of  200  -  300,  or  even  higher,  depending  on 
other  physics  constraints2.  This  suggests  that  high- aspect  ratio  shells  are  indeed  a 
viable  and  highly  attractive  design  for  direct-drive  laser  fusion.  Clearly,  we  now  need 
additional  experimental  data  on  the  RT  instability  under  various  conditions  to  test 
these  predictions.  We  thank  P.  C.  Reed  for  her  technical  assistance.  This  work  was 
supported  by  the  U.  S.  Department  of  Energy  and  the  Office  of  Naval  Research. 
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Figure  1.  Comparison  of  the  classical  growth  rate  ((kg)1'*, — ■)  to  the  simulation  results 
(•  )  and  to  the  theoretical  result  (Eq.  1,  x)  for  1/im  laser  light  with  an  absorbed  intensity 
of  1013  W/cm2  impinging  20^m  thick  plastic  (CH)  foils.  The  ablation  velocity  is 
measured  at  the  center  of  the  shear  layer,  O(lum)  from  the  target  surface,  and  is  the 
same  for  ail  perturbation  wavelengths,  \\  =  1.3  x  I0°cm,  s  . 


Figure  4.  Isovorticity  contours,  in  10  %  increments  of  the  maximum,  on  the  target 
surface  eariy  in  the  evolution  of  the  20^m  RT  instability  for  the  l^m  laser  case.  The 
laser  is  impinging  the  target  from  the  right.  The  spatial  dimensions  are  in  fj,i n.  The 
dashed  (solid)  contours  indicate  clockwise  (counter  clockwise)  flow.  The  wavy  bold 
lines  deiineate  the  10  %  (right)  and  SO  %  (left)  density  contours.  The  ablation  velocity 
(Va,  em.'s),  which  is  convecting  vorticity  away  from  the  target  surface,  is  plotted  on  top 
of  the  vorticity  contours.  The  value  of  at  the  center  of  the  layer  is  1.3  x  I0°cm/s. 


Figure  5.  Plot  of  the  ablation  velocity  (Vi,)  as  a  function  of  distance  from  the  target  sur¬ 
face  (Y)  for  1/zra  (solid),  kfim  (dash)  and  (dash-dot)  laser  light  at  3  x  10 IZW ;  cm? 
impinging  a  flat  20^m  thick  piastic  (CH)  foil.  The  curves  begin  to  cross  at  distances 
0(40um)  from  the  target  surface.  The  vertical  lines  bracket  the  vortex  layer  (see  Ref. 


References 


1.  Lord  Rayleigh,  Theory  of  Sound  (Dover  Publications,  Inc.,  New  York  (1894),  2nd 

ed.,  Vol.  2;  G.I.  Taylor,  Proc.  Roy.  Soc.  (London)  A201.  192  (1950). 

2.  S.E.  Bodner,  J.  Fus.  Energy  i  221  (1981). 

3.  M.H.  Emery  et  al.,  Appl.  Phys.  Lett.  41,  808  (1982). 

4.  J.H.  Gardner  and  S.E.  Bodner,  Bull.  Am.  Phys.  Soc.  29,  1290  (1984);  S.E.  Bodner 

et  al..  Plasma  Phys.  and  Cont.  Nuc.  Fus.  Res.  3,  155  (1984). 

5.  C.P.  Verdon  et  al.,  Phys.  Fluids  25,  1653  (1982);  M.H.  Emery  et  al.,  Phys.  Rev. 
Lett.  48,  677  (1982). 

6.  R.R.  Whitlock  et  al.,  Phys.  Rev.  Lett.  52,  819  (1984);  J.  Grun  et  al.,  Phys.  Rev. 

Lett.  53,  1352  (1984). 

7.  M.H.  Emery  et  al.  Phys.  Rev.  Lett.  48,  253  (1982);  P.G.  Burkhalter  et  al.,  Phys. 

Fluids  26,  3650  (1983);  M.J.  Herbst  et  al.,  Phys.  Rev.  Lett.  52,  192  (1984). 

8.  M.H.  Emery  et  al.,  Phys.  Fluids  27,  1338  (1984). 

9.  M.H.  Emery,  Bull.  Am.  Phys.  Soc.,  Invited  paper  3A5,  29,  1231  (1984);  to  be 
submitted  to  Phys.  Fluids. 

10.  W.M.  Manheimer  and  D.G.  Colombant,  Phys.  Fluids  27,  983  (1984). 

11.  J.H.  Gardner  and  S.E.  Bodner  Phys.  Rev.  Lett.  47,  1137  (1981). 


